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Abstract. 

Phase equilibria between regions of different thickness in thin liquid films stabilized by 
colloidal particles are investigated using a quasi-two-dimensional thermodynamic formalism. 
Appropriate equilibrium conditions for the film tension, normal pressure, and chemical potential 
of the particles in the film are formulated, and it is shown that the relaxation of these parameters 
occurs consecutively on three distinct time scales. Film stratification is described quantitatively 
for a hard-sphere suspension using a Monte-Carlo method to evaluate thermodynamic equations 
of state. Coexisting phases are determined for systems in constrained- and full-equilibrium 
states that correspond to different stages of film relaxation. 



Drainage of thin liquid films often involves spontaneous formation of coexisting regions 
of different, but uniform, thickness. For example, circular black-film spots can be seen in 
soap bubbles prior to breakup. Such a stratification phenomenon is particularly interesting in 
films stabilized by colloidal particles, micelles, or polyelectrolytes, where the thinning occurs 
through a series of stepwise transitions between film states characterized by the thickness 
commensurate with the size of the stabilizing particles. In experiments with horizontal films 
[1,2], such transitions occur by rapid formation of circular regions of a smaller thickness, 
followed by a much slower expansion. In investigations of vertical films, up to seven coexisting 
parallel stripes of different thickness have been observed [3]. 

Film stratification has recently been intensively studied [1-8] because of the relevance of 
the problem for understanding structural colloidal forces. In particle-stabilized thin liquid 
films, the structural forces depend on the film thickness in an oscillatory manner due to par- 
ticle layering [9] . Equilibrium values of the film thickness result from the balance between the 
normal pressure in the film (which includes the oscillatory structural force) and the outside 
pressure [6,7,10]. We find, however, that the normal-stress balance used alone in descriptions 
of the film stratification phenomenon is insufficient — the mechanical equilibrium between re- 
gions (phases) of different thickness requires also the lateral force balance. Since the lateral 
balance has not been included in the available theories, understanding of an essential aspect of 
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the problem is still lacking. One of our goals is to provide a proper thermodynamic framework 
for analyzing the phase behavior of particle stabilized thin liquid films. 

An analysis of the stepwise-thinning phenomenon requires equilibrium as well as nonequi- 
librium considerations. Experimental observations indicate that the thinning process occurs 
on several distinct time scales: a thinner spot nucleates in a fraction of a second but the time 
scale for the subsequent expansion of the spot is much slower. This separation of the time 
scales suggests that the system at different evolution stages is, approximately, in constrained 
equilibrium. In such states the film phases of uniform thickness arc in equilibrium. However, 
only some of the equilibrium conditions between these phases are satisfied — the others are 
not, due to dynamical constraints. An important goal of our study is to identify relevant 
time scales in the stepwise-thinning phenomenon, and to discern the corresponding relaxation 
processes. Our theoretical considerations are supplemented with results of Monte-Carlo simu- 
lations of a model system to illustrate possible full and partial equilibrium states of a stratified 
particle-stabilized film at different evolution stages. 

Our analysis of the film stratification phenomenon relies on representing the film as an 
effective two-dimensional medium [9, 11]. Accordingly, all details of the transverse structure 
of the film are averaged out (as in the thermodynamic description of an interface proposed 
by Gibbs). The film thickness, however, is retained as a parameter of state, because it can be 
controlled by varying the difference between the external pressure p^ and the pressure of the 
bulk fluid in contact with the fluid in the film. We note that the usual formulation of thin-film 
thermodynamics [9] is insufficient for our problem because of the underlying assumptions that 
the pressure in the film is isotropic and the interaction between the film interfaces is described 
by a mean-field normal force. In contrast, the stabilizing colloidal particles contribute not 
only to the normal but also to the tangential stress balance in the film, and the corresponding 
stress tensor is non-isotropic due to confinement effects. These features have to be taken into 
account in the proper formulation of the theory. 

A typical thickness of a particle-stabilized liquid film is in the micrometer range, consis- 
tent with the size of colloidal particles. On this scale, the film interfaces can be treated as 
surfaces with zero thickness and interfacial tension a. Assuming this geometry, a suitable 
effective two-dimensional thermodynamic description can be obtained from the expression for 
the macroscopic infinitesimal work 



associated with changes of the film thickness h and area A. In the above relation, p± and pu 
denote the normal and lateral components of the pressure tensor in the film and the overbar 
indicates the average value across the film. If a direct interaction force (e.g., the van der Waals 
attraction or electrostatic repulsion) acts between the film interfaces, the normal pressure p± 
in eq. Q should be replaced with p± + /„, where f n is the normal force per unit area, and 
/„ > corresponds to attraction. 

In eq. Q , the work associated with the area change consists of the interfacial-tension and 
pressure contributions. To obtain an effective two-dimensional thermodynamic description 
of the film, the lateral interfacial-tension forces and the pressure are combined into a single 
conjugate thermodynamic variable. However, first an appropriate set of independent thermo- 
dynamic parameters has to be chosen. The pair of independent variables h and A used in Eq. 
fl} is unsuitable because h is not an extensive quantity. By a simple change of variables, a 
more appropriate form is obtained 



where V = Ah is the film volume, and, by analogy with the definition of the surface tension, 



dW = -Apx dh + (2a - hp\\ ) dA 



(1) 



aw = -p± dV + jdA, 



(2) 



7 = h(p± ~P\\) + 2cr 



(3) 
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is the tension of the film [11]. 

The fundamental relation for the film free energy F can be derived in the usual way from 
the expression for mechanical work J2J • Treating the suspension in the film as a two-component 
fluid (with the colloidal particles regarded as macromolecules) we get 

dF = -SdT- Pl _ dV + 7dA, +/i c dN c + Mf dN { , (4) 

where N c and Nf denote the number of colloidal particles and solvent molecules in the system, 
and fi c and /if are the corresponding chemical potentials. The film-tension representation 
Q of the film free energy is more suitable for discussion of phase equilibria in stratified 
films than the more common film surface tension representation where the film thickness h is 
retained as an independent thermodynamic parameter. Since all independent parameters in 
our representation V, A, N c and Nf are extensive, the conjugate intensive functions of state 
p^,7,/i c , and /if in coexisting phases must be equal in unconstrained equilibrium (which 
follows from standard entropy-maximization arguments [12]). Accordingly, apart from the 
trivial thermal-equilibrium condition = = T^ c \ there are the transverse and lateral 
mechanical equilibrium conditions 

p W =p ( L 2 ) =p (e), 7 (1) =7 (2), (5a , b) 

and the chemical equilibrium conditions 

/4 X) =/4 2) , /4 1} =M f 2) , (6a,b) 

where and are the external temperature and pressure. 

A typical size of colloidal particles used in experiments with stratified thin films (or the 
range of the effective screened Coulombic potential if the particles are charged) is > 0.1 /im. 
Under these conditions the solvent can be treated as a structureless incompressible continuous 
medium that does not affect thermodynamic properties of the system. In our further consid- 
erations the solvent variables /if and Nf are thus neglected (which is a usual approximation), 
and the equilibrium conditions (J5J arc rephrased in terms of the excess quantities associated 
with the particle presence, 

P2=P2, 7«=7< 2) , (7a,b) 

where the index c indicates the particle contributions. 

As demonstrated in numerous experiments, formation of a stepwise structure in particle- 
stabilized films is much faster than the subsequent evolution of the coexisting phases [1-3]. 
Scaling arguments, outlined below, suggest that the relaxation of a stratified film is charac- 
terized by three distinct time scales tii <C ti < t c for attaining the equilibrium conditions 
Cb)j Ct)i an d respectively. The corresponding noncquilibrium thermodynamic pro- 

cesses entail changes in the areas of individual phases, their volumes, and the volume fraction 
of colloidal particles. 

Estimates of the characteristic relaxation times can be obtained using energy-dissipation 
arguments. To evaluate the time scale tii we note that the area change involves energy 
dissipation in the whole volume of the film V and is characterized by the velocity variation on 
the lengthscale I set by the lateral film dimension. By comparing the energy-dissipation rate 
Vrj{u/l) 2 to the power ulhn c kT associated with the thermal stresses n c kT produced by the 
particles, we find tii ~ rj/n c kT, where ij is the fluid viscosity, u ~ l/rn is the characteristic film 
velocity, n c = N c /V, and kT is the thermal energy. The change in the volume of the coexisting 
phases requires suspension flow through the contact region where the film thickness changes 
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Fig. 1 - Particle contribution to normal pressure p c ± and film tension y c versus film thickness h for 
particle volume fraction 4> c = 0.4. Dotted lines correspond to states for which the stability condition 
HUH is violated. Circles indicate a pair of stable phases in mechanical equilibrium Q. 



rapidly. Accordingly, the magnitude of the velocity gradient is 0{u/h) and the dissipation 
occurs in the volume (h/l)V, which yields the estimate t± ~ (l/h)n . A change of the particle 
volume fraction </> c involves particle motion with respect to the fluid in the whole volume of 
the film. Assuming d ~ h (where d is the particle diameter) and 4> c — 0(1) we find that 
t c ~ (l/h)r±. For an aqueous suspension of particles with the diameter d = 0.1 /jm in a film 
with the characteristic lateral dimension / = 10 /im we find tii ~ 10~ 3 s, tj_ ~ 10 _1 s, and 
t c ~ 10 s, which is consistent with experimental observations [2]. We note that our scaling 
analysis applies to films with surfactant-free interfaces (as in the experiments reported in [2]). 
We expect that a surfactant adsorbed on film interfaces may increase the time scale t± by a 
factor l/h, due to the modification of hydrodynamic boundary conditions. Also, the timescale 
t c may be decreased by a factor d/h for charged particles with d <C fa- 
Assuming the separation of time scales tii < ti C t c we anticipate the following scenario 
for the stepwise-thinning of the film. At an early stage of phase separation, the relaxation 
occurs on the time scale t\\ primarily through a change in the areas of the coexisting phases at 
constant <j> c . The process is driven by an imbalance of lateral forces, and after it is completed 
the film is in a constrained equilibrium state where only the lateral-force equilibrium condition 
|0d ) between the coexisting phases is satisfied. (An analogy is the coexistence of water and 
air where the pressure difference between the two phases relaxes very fast, but equilibration 
of the temperature and chemical-potential differences is much slower.) The next stage of the 
film evolution is driven by the difference Ap c ± of the colloidal contribution to the normal 
pressure in the coexisting regions. As required by the incompressibility of the suspension in 
the film, this pressure difference is compensated by the dynamic pressure drop due to the 
suspension flow through the contact region. After the pressure difference Ap c ± has relaxed on 
the time scale t± the chemical potential of colloidal particles relaxes on the time-scale t c due 
to the particle diffusion in the film. Only this last process entails a change of particle volume 
fraction in the coexisting regions. 

Our theory can be used to obtain quantitative predictions for the trajectory in the ther- 
modynamic parameter space for a stratified film during the equilibration process. To this end 
the general thermodynamic relations have to be supplemented with equilibrium equations of 
state. We present numerical results for a film stabilized by a monodisperse suspension of 
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Fig. 2 - Chemical potential of the particles fi c versus film thickness h for <j> c = 0.4. The normalization 
corresponds to the asymptotic behavior fi c — — \nv for 4> c <C 1. Dotted line as in fig.Q 

hard spheres. The particle contribution to average stress tensor p~ c in homogeneous phases 
is evaluated from the contact value of the equilibrium two-particle density n^ri,^) in the 
confined suspension using the relation 

£=<+ (/ frn^( ri ,n - dr) d 2 f^ . (8) 

Here f = r/r, the integration is over the contact configurations, (• • • )y denotes the volume 
average over the region V = A(h — d) accessible to particle centers, and n' c — N c /V is the 
corresponding particle number density( 1 ). Equation (JSJ generalizes a well-known expression 
for the pressure in a bulk hard-sphere system in terms of the contact value of the radial 
distribution. It can be derived either from the collisional contribution to the momentum flux 
or by passing to the hard-sphere limit of the Kirkwood-Buff expression [13] for the stress 
tensor in an inhomogeneous fluid. The equilibrium average involved in JHJ was evaluated 
via computer simulations using the standard Metropolis algorithm to generate the canonical 
ensemble for a system of several hundred spheres confined in a planar film with periodic 
boundary conditions in the lateral directions. 

A typical dependence of the normal pressure p c ± and film tension j c on the film thickness h 
at a constant volume fraction </> c is shown in fig.^ By integrating the Gibbs-Duhem relation 
at fixed T 

d^ c = v dp c± - a djc (9) 

(where v — V/N c and a = A/N c ), we also determined the chemical potential /i c . The plot of 
/i c versus h at constant <p c is presented in fig.|21 The film thickness h is normalized here by the 
particle diameter d, pressure p c ± by kT/d 3 , film tension j c by kT/d 2 , and chemical potential 
He by kT. As expected, the intensive parameters p c ±, 7 C , and /i c oscillate as functions of the 
film thickness h. Simulation results at different volume fractions (not shown) indicate that 
the oscillations are more pronounced at higher <j) c . We note that the thermodynamic stability 
of the system requires that 

(d Pc Jdh) Tllc < 0, (10) 

( 1 )The accessible volume V' rather than the total volume V appears in eq. J^J because p" c is the average of 
the local stress tensor p c (i") that vanishes for r V' . Consistently, we use the relation 7 C = h'(p~^ — Pc±) for 
the particle contribution to the film tension. 
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Fig. 3 Fig. 4 

Fig. 3 - Parametric plot of the equation of state in the (p c x,7c) plane at constant particle volume 
fraction <j) c — 0.4. The curve is parametrized by film thickness h; direction of the increasing h is 
indicated by the arrow. Dotted line as in fig. Q Circle indicates the stable coexistence state. 

Fig. 4 - Same as fig. 01 except that the results are for the constant chemical potential fi c — 11. Inset 
shows the variation of the particle volume fraction with the film thickness h along the curve. Circles 
indicate stable coexistence states; film thickness in the coexisting phases is (/i.1,/12) — (1-1,3.0), 
(1.1,2.1), and (2.0,3.0). 



which is equivalent to (dj c /dh)T^ c < 0. The domains where Hl()[) is violated are indicated in 
the figures by dotted lines. 

We now return to our discussion of phase equilibria at different stages of film evolution. 
As predicted by our scaling analysis, the first two stages occur at a constant volume fraction 
4>c^ = <pc 2 ^ ■ At the end of the first stage, the coexisting phases fulfill the lateral mechanical- 
equilibrium condition (7p); accordingly, the film thickness in these phases satisfies equation 
7 c (/i) = 70, where the film tension is evaluated at constant (ft c , as in fig.^. The value of the 
constant 70 is not determined by the equilibrium condition (J2)) alone, but it also depends on 
the initial volumes of the coexisting regions. The oscillatory character of the curve 7 C = 7c (h) 
in fig. indicates that several different combinations of coexisting phases are possible. 

At the next relaxation stage, the complete mechanical equilibrium is attained. Since two 
equilibrium conditions J7J) are now satisfied, the coexisting phases correspond to isolated points 
in the space of states for given </> c . The coexisting phases can be determined from the equations 
of state p c ±_ = p c j_(h) and 7 C = ^ c {h) by combining them into a single parametric plot, as 
shown in fig.|3for (j> c = 0.4. The coexisting phases correspond to the line intersections. In our 
example only one intersection (indicated by the circle) represents a stable system. The two 
coexisting phases are also marked in fig. Additional numerical calculations show that at 
higher volume fractions several distinct stable states with coexisting phases can be identified, 
and for lower <p c there are no such states. We emphasize that not only do the intersection points 
in fig.|21have a physical meaning but also the curve itself does — its different portions represent 
the trajectory of the system in the space of states when the normal-pressure difference relaxes. 

A similar parametric plot can be used to determine the possible coexisting phases in full 
thermodynamic equilibrium. However, the chemical potential /i c is now held constant, instead 
of </> c , due to the equilibrium condition JSJi). An example of such a plot is presented in fig. 01 
the inset shows the corresponding variation of the volume fraction </> c with the film thickness 
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h. The shape of the curve /i c = const in the (p c _L,7c) plane is entirely different from the 
shape of the line 4> c — const in fig. [31 The distinct qualitative features of the curve plotted 
in fig. 0] follow from Frumkin [14] equation (dj c /dp c ±)Tp, c — h (which can be derived from 
(jHJ). Frumkin equation implies that the slope of the line j c = 7c(pcJ_) at constant /i c is always 
positive and grows with the increasing h. Moreover, the extrema of the functions p c ±(h) and 
j c (h) at constant /x c coincide, which causes the sharp corners in the plot. 

The line j c = 7 c (p c _l) in fig-E]has numerous intersections, indicating a number of possible 
combinations of coexisting phases. In contrast, the results in fig.[3]have only one intersection. 
For a given system several equilibrium states may thus exist, but it may be difficult to reach 
them in a stepwise-drainage process because of a bottleneck at the second relaxation stage: 
the evolving coexisting phases may encounter an instability before mechanical equilibrium (JJJ 
is reached. After one or more of such transitions the film may break rather than achieve a 
stable unconstrained-equilibrium configuration. (Indeed, film breakup was observed in most 
experiments with particle-stabilized films; see, however, recent results reported in [2].) 

Designing experiments that would clearly distinguish between different evolution phases 
of the stepwise-drainage process is a challenge that requires further numerical and theoretical 
studies. Our approach will be used in our forthcoming study to propose experimental protocols 
for reaching full equilibrium states predicted by our analysis. 

* * * 
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